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Abstract. We show that the multiplication operator associated to a fractional power of a 
Gamma random variable, with parameter q > 0, maps the convex cone of the 1-invariant 
functions for a self-similar semigroup into the convex cone of the q-invariant functions for the 
associated Ornstein-Uhlenbeck (for short OU) semigroup. We also describe the harmonic func- 
tions for some other generalizations of the OU semigroup. Among the various applications, we 
characterize, through their Laplace transforms, the laws of first passage times above and over- 
shoot for certain two-sided a-stable OU processes and also for spectrally negative semi-stable 
OU processes. These Laplace transforms are expressed in terms of a new family of power series 
which includes the generalized Mittag-Leffler functions and generalizes the family of functions 
introduced by Patie [20f | - 



1. Introduction and main results 

Let E = R, R + or [0, oo) and let X be the realization of (Pt)t>o, a Feller semigroup on E 
satisfying, for a > 0, the a-self-similarity property, i.e. for any c > and every / G 2J(E), the 
space of bounded Borelian functions on E, we have the following identity 

(1) P d * t f(cx)=P t (d c f)(x), xGE, 

where d c is the dilatation operator, i.e. d c f(x) = f(cx). We denote by (F x ) x& e the family of 
probability measures of X which act on D(E), the Skorohod space of cadlag functions from 
[0, oo) to E, and by (F*)t>o its natural filtration. We also mention that, throughout the 
paper, E stands for a reference expectation operator. Moreover, A (resp. 2) (A)) stands for its 
infinitesimal generator (resp. its domain). We have in mind the following situations 

(1) E = R and X is an a-stable Levy process. 

(2) E = R + or [0, oo) and X is a —semi-stable processes in the terminology of Lamperti. 

More precisely, let £ be a Levy process starting from x S R, Lamperti [l6|] showed that the time 
change process 

(2) X t = e^t , t > 0, 
where 

A t = mt{u > 0; V u := / ds > t}, 

Jo 
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is an —semi-stable positive Markov process starting from e x . Further, we assume that — oo < 
b := E[£i] < oo and denote the characteristic exponent of £ by ^. We also suppose that £ is 
not arithmetic (i.e. does not live on a discrete subgroup rTL for some r > 0). Then, on the one 
hand, for b > (resp. b > and £ is spectrally negative), it is plain that X has infinite lifetime 
and we know from Bertoin and Yor 0, Theorem 1] (resp. 0, Proposition 1]), that the family 
of probability measures (¥ x ) x> q converges in the sense of finite-dimensional distributions to a 
probability measure, denoted by Po, as x —* 0+. On the other hand, if b < 0, then it is plain that 
£ drift towards — oo and X has a finite lifetime which is = inf{n > 0; X n _ = 0, X u = 0}. In 
this case, we assume that there exists a unique 9 > which yields the so-called Cramer condition 

(3) E[e e?1 ] = 1. 

Then, under the additional condition < 9 < a, Rivero 24] showed that the minimal process 
(X,T X ) admits an unique recurrent extension that hits and leaves continuously a.s. and which 
is a —semi-stable process on [0, oo). With a slight abuse of notation, we write (f x )x>o for the 
family of laws of such a recurrent extension. We gather the different possibilities in the following. 



HO. E = [0, oo) and either b > or b < and ((3} holds with 9 < a. Moreover, if £ is 
spectrally negative, the case b = is also allowed. 



We will also use the following hypothesis 



HI. E C [0, oo). £ has finite exponential moments of arbitrary positive orders, i.e. ip(jn) < oo 
for every m > where ip(m) = —^(—im). 



If £ is spectrally negative, excluding the degenerate cases, then HI holds and for b £ R, we 
write 9q for the largest root of the equation t/j(u) = 0. Then, being continuous and increasing 
on [#o,oo), tp has a well-defined inverse function <f> : [0, oo) — > [9q,oo) which is also continuous 
and increasing. 

For r > 0, we write P[ = e~ rt Pt- We say that a non-negative function, I r , is r-excessive 
(resp. r-invariant) for P± if for any t > 0, 

P[l r (x) < J r (x), 

(resp. if we have = in place of <) and limtio P[I r (x) = 2 r (x) pointwise. Taking r = 1 and 
writing simply X = T\, we have P}l(x) = T{x). The self-similarity property ([1]) then yields 

(4) /f(dil)(r4x)=J(i), 

which entails the identity Z r {x) = 2(rax) for all x £ E and r > 0. We denote the convex cone 
of 1-excessive (resp. 1-invariant) functions for X by <£(X) (resp. 3(X)). 

For any A > 0, the Ornstein-Uhlenbeck (for short OU) semigroup, (Qt)t>o, is defined, for 
/ e 55(E), by 

(5) Qtf(x) = P ex{t) (d e ' x (-t)f) (x), xGE,t>0, 

where e\(t) = e V" 1 , x = a ^ an d we write v\(.) for the continuous increasing inverse function of 
e\(.). We mention that such a deterministic transformation of self-similar processes traces back 
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to Doob who studied the generalized OU processes driven by symmetric stable Levy pro- 
cesses. Moreover, Carmona et al. 0, Proposition 5.8] showed that (Qt)t>o is a Feller semigroup 
with infinitesimal generator, for / G D (A), given by 

U/0) = A/(x) - Xxf'(x). 

Let U be the realization of the Feller semigroup (Qt)t>o- It follows from ([5]) that 

(6) U t = e' x (-t)X eAt) ,t>0. 

We deduce, with obvious notation, that Tq = v x (T^) a.s.. If E = M (resp. otherwise), we call 
U a self-similar (resp. semi-stable) OU process. We denote by (Q x )x>o the family of probability 
measures of a semi-stable OU process. We deduce, from the Lamperti mapping ([2]) and the 
discussions above the following. 

Proposition 1.1. For any x > 0, there exists a one to one mapping between the law of a Levy 
process starting from log(x) and the law of a semi-stable OU process starting from x. More 
precisely, we have 

(7) U t = e~ xt e^, t <T C/ , 

where At = Jjj U~ a ds. Note that for b > 0, the previous identity holds for any t > 0. 

Moreover, if ([3]) holds with < 6 < a, then the minimal process (U,Tq) admits a recurrent 
extension which hits and leaves continuously a.s. which is the OU process associated to (X,¥ x ). 
We write its family of laws by (Q x ) x >o- 

Under the condition HI, for any 7 > 0, we denote by the law of the semi-stable OU 
process, starting at x G M + , associated to the Levy process having Laplace exponent ^ 7 (u) = 
ip(u + 7) — ^(7), u > 0. 

We say that a probability measure m is invariant for U if it satisfies, for any / G 53(E), 



Qtf(x)m(dx) = / f{x)m{dx). 
Je 

Finally, let G q be a gamma random variable independent of X, with parameter q > 0, whose 
law is given by j(dr) = e dr. We are now ready to state the following. 

Theorem 1.2. If E = ~R or HO holds then the Feller process U is positively recurrent and its 
unique invariant measure is x^o(Xi G dx). 

Next, assume that I G L 1 (7(dr)). For any q > 0, we introduce the function I (q; x) defined by 



(8) J(g;x)= X *E 



1 (XG 



x G E. 



Then, ifl G 3(X) (resp. <£(X)) then T{q;x) G 3 q (U) (resp. £ q {U)). 
Consequently, ifZ£ 3(X), we have, for any q > 0, 

(9) (l+xtriPt(d (i+xtr xl)(q;x)=l(q;x), xeE. 

Remark 1.3. (1) We call the multiplication operator ([8]) associated to a fractional power of 
a Gamma random variable, the T-transform. 
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(2) The characterization of time-space invariant functions of the form Q, associated to 
self-similar processes, has been first identified by Shepp [2(i ] in the case of the Brown- 
ian motion and by several authors for some specific processes: Yor 23] for the Bessel 
processes, Novikov fl9[ ] and Patie 22] for the one sided-stable processes. Whilst in the 
mentioned papers, the authors made used of specific properties of the studied processes 
to derive the time-space martingales, we provide a proof which is based simply on the 
self-similarity property. 

We proceed by investigating the process Y, defined, for any x,/3 G R and Co = a.s., by 
(10) Y t = e a ^ (x + P f e'^ds] ,t>0. 



We call Y the Levy OU process. We mention that this generalization of the OU process is a 
specific instance of the continuous analogue of random recurrence equations, as shown by de 
Haan and Karandikar They have been also well-studied by Carmona et al. [fj], Erickson and 
Mailer [13] , Bertoin et al. [l[ and by Kondo et al. 14 ] . In 0] , it is proved that Y is a homogeneous 
Markov process with respect to the filtration generated by £. Moreover, they showed, from the 
stationarity and the independency of the increments of £, that, for any fixed t > 0, 



Y t = xe a ^ + P f e^'ds. 
Jo 



Then, if E[£i] < 0, they deduced that, as t — ► oo, & —>■ — oo and Y t — > PVoo = Jq°° e a ^ s ds. We 
refer to Bertoin and Yor [4J for a thorough survey on the exponential functional of Levy processes. 
In the spectrally negative case, it is well know that the law of is self-decomposable, hence 
absolutely continuous and unimodal. Moreover, under the additional assumption that 9 < a, 
its law has been computed in term of the Laplace transform by Patie [20]. Now, we introduce 
the process Z defined, for any x ^ 0, [5 £ R and £o = a.s., by 

(11) Z t = e a ^ (x + P f e^ds] , t > 0. 



Before stating the next result, we introduce some notation. Let B be a Borel subset of E and 
we write Tg for the first exit time from B by U . With a slight abuse of terminology, we say 
that for any x 6 E, a non-negative function H is a (qA, S)-harmonic function for (U, <Q X ) if 

(12) K x [e~ qATl B-H{U T u)I {T v <T u } ] = H(x). 

When At is replaced by t in the previous expression, we simply say that Ti is a (q, 5)-harmonic 
function for (U,Q X ). We are ready to state the following. 

Theorem 1.4. Set P = aXx in (|11|) . Then, to a process Z starting from ^ with i / 0, one 
can associate a semi-stable OU process (U,Qi) such that 

(13) Z t = x~ 1 U^ t ,t<T^, 

where \/t = f Z s ds and its inverse is given by At = Jq U^~ a ds. Note that for b > 0, the previous 
identity holds for any t > 0. Consequently, with x > 0, Z is a Feller process on (0, oo). 

Moreover, let q > ; < a < b < +oo and x G (a,b). Then, a (q A, TF 6a .p" harmonic function 
for (U, Qi) is a (q, T? a <x \,<x\) -harmonic function for the process Z starting from x~ a . Similarly, a 



(qA,TK x\) -harmonic for (U, Qi) is a (q,TX a , a ^) -harmonic function for the process Y starting 

from x a , the Levy OU process associated to the Levy process £ = — £, the dual of £ with respect 
to the Lebesgue measure. 

Finally, assume that HI holds and write p q {x) = x g for x,q > 0. // the function 7i is 

(Xcj)(q), B) -harmonic function for (E/,Qas ) then the function p^^TL is (qA, B)-harmonic func- 
tion for (U,Q X ). 



2. Proofs 



2.1. Proof of Theorem 11.21 The description of the unique invariant measure is a refinement 
of 0, Proposition 5.7] where therein the proof is provided for R-valued self-similar processes 
and can be extended readily for the K + - valued case under the condition HO, which ensures that 
(X, ¥ x ) admits an entrance law at 0. 

Next, let us assume that X G L 1 { r ^(dr)) n 3(X). We need to show that for any q > 0, 
e- qt Q t l(q;x) =l(q;x). For iEE, we deduce from the definition of (Qt)t>o that 



e- gt Q t l(q-x) 



X x 



X x 



J x((x r ) a Utje r rx l dr 
l((xr)«e> x (-t)X ex{t) ) 



e r x dr 



Using the change of variable u = X e v( — *) r i Fubini theorem and ©, we get 



e- qt Q t l{q ]X ) 



E, 



| [£) ./n 



e x (*) 



e xux du 



I ( u a x ) e *ti* 1 c?u 



= 2% af) 

where the last line follows after the change of variable u = x r - The case I € L 1 (7((ir)) n £(X) 
is obtained by following the same line of reasoning. The last assertion is deduced from ([5]) and 
(jHJ) by performing the change of variable u = v x (t), with v x (t) = ^ log(l + x~k)- 



2.2. Proof of Theorem 11.41 Setting (3 = aXx, the Lamperti mapping 

Z t = x~ x e a ^ ( 1 + aX / e a ^ds 



yields 



o 



1 



(1 + aX.) 



-x. 



a 



V t 



x- 1 (a 



v x(-))v t 
5 



where the last identity follows from © . The proof of the assertion (|13p is completed by observing 
that 

-t \ -i 



(v x (V t ))' = e a ^ (l + axj e a ^ds 



Moreover, since the mapping x i— > x a is a homeomorphism of M + , the Feller property follows 
from its invariance by "nice" time change of Feller processes, see Lamperti fl~5l . Theorem 1]. We 
also obtain the following identities 

T ( f )6) = inf{u>0; Z u ${a,b)} 

= inf {u > 0; U* u i {ax,bx)} 

= A(inf{n>0; U u $ ((aa;)£, (6z)£)}) . 

The characterization of the harmonic functions of Z follows. The characterization of the har- 
monic functions of Y are readily deduced from the ones of Z and the identity 

Y t =^r,t>Q. 

The proof of Theorem is then completed by using the following Lemma together with an appli- 
cation of the optional stopping theorem. 

Lemma 2.1. Assume that HI holds, then for 7, 5 > and x > 0, we have 



(1) d Q(7) = ^V 5 e X(y-S)t- Wl )- m )A tdQ (8)^ onF U n { t<T V}. 

Note that for b > the condition " on {t < Tq}" can be omitted. For the particular case 7 = 6 
and 5 = 0, the absolute continuity relationship (pQ) reduces to 

w = (EiY e ™t dqxj onF u n{t<T u } 



t.r 



X 



Proof. We start by recalling that in [20(, the following power Girsanov transform has been 
derived, under HI, for 7, 5 > and x > 0, with obvious notation, 

eflpW = ^y~Vw(7)-^))^ rfp(f), on F t x n {t < }. 

The assertion ([!]) follows readily by time change and recalling that 



A ex (t) = / U~ a du. 

Jo 

We complete the proof by recalling that for b > 0, U does not reach a.s.. □ 
3. Applications 

In this section, we illustrate our results to some new interesting examples. 
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3.1. First passage times and overshoot of stable OU processes. Let X be an a-stable 
Levy process whose characteristic exponent satisfy, for m£|, 



Vl/(wt) = — c\u\ a ^1 — i/3sgn(u) tan 



where 1 < a < 2 and for convenience we take c = (1 + (3 2 t&n 2 (jf)) 1/2 . Then, we introduce 
the constant p = P(Xi > 0) which was evaluated by Zolotarev [28| as follows 



P 



- + — tan -1 ( P tan ( ) 



na 



V 2 J J 



Following Doney 10], we introduce, for any integers k,l, the class C& j of stable processes such 
that 

p + k = la 

where a = ^. For m € N, x G M and z£C, introduce the function 



f m (x,z) = Y[ (z + e 



J,x(m—2i)TT 



i=Q 



Next, we recall from the Wiener- Hopf factorization of Levy processes due to Rogozin [25|, that 
the law of the first passage times and the over(under)shoot of X at the level is described 
by the following identities, for 5, r > and p > 0, 



f 

Jo 



e~ 5x E. 



e~ Sx E x 



-rr^-pX x _ 
e 



-rr^-pX x 

e o 



dx 



dx 



1 



5 — p 
1 

5 — p 



$+(-r«tf) 
$+(-rop) / 

^-(r«p) J 



where (1 - ^O*)) -1 = tf-(<5)tf+(<5). Here *+(<5) (resp. #~(<5)) is analytic in < (resp. > 
0) continuous and nonvanishing on "Si(5) < (resp. ?ft(5) > 0). Doney [10j] computes the Wiener- 
Hopf factors for stable processes in i as follows 



* + (z) 



/ w (5,(-l)^+^) 

/;_!(«, (-l) fc +^) 

/ fc (a, (-!)'*«) ' 



, atfl(z) / o, 

2ltg(z) / -7r, 



where z^ 3 stands for a^e 1 ^ when z = ae 1 ^ with u > and — tt < (b < it. Observe also that 



^ + (— xa ) ~ x p for large real x. Moreover, using the fact that the function E a 



e o 



is r-excessive for the semigroup of X, we deduce from the T-transform the following. 
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Corollary 3.1. For any q, 5 > 0, p > 0, and for any integers k,l such that X £ C^i, we have 



o 



e Sx E* 



oo 
oo 



e o 



e o 



dx 



dx 



1 



5 — p 
1 

5 — p 



X x 



1 



r (f ) Jo v+(-Ap) 
i 



e x r x 



r (f ) ■/ o ^~(r«p) 



e x r x ctr 



3.2. First passage times of one-sided semi-stable- and Levy-OU processes. We now 

fix (Pt)t>o to be the semigroup of a spectrally negative —semi-stable process X. X is then 
associated via the Lamperti mapping |2| to a spectrally negative Levy process, £, which we 
assume to have a finite mean b. Its characteristic exponent if) has the well known Levy-Khintchine 
representation 



(1) 



ip(u) =bu + —v 2 + / 

^ J-oc 



(e ur - \-ur)v{dr), u > 0, 



where a > and the measure v satisfies the integrability condition J ^^ A r 2 ) u(dr) < +oo. 
Patie 20] computes the Laplace transform of the first passage times above of X as follows. For 
any r > and < x < a, we have 



(2) Ex ■ T , 

where the entire function, I a ^, is given, for 7 > and a > 0, by 



n=0 



and 



a n (ip;a) 1 = ip(ak), a = 1. 
fc=i 

Using the T-transform, we introduce the following power series 

00 

(3) 2a,V>(<?; Z) = ffln ^' a )(^)n^ n 

n=0 

where ((;) n = is the Pochhammer symbol and we have used the integral representation 

of the gamma function T(q) = e~ r r q ~ 1 dr, $i(q) > 0. By means of the following asymptotic 
formula of ratio of gamma functions, see e.g. Lebedev [13, p. 15], for 5 > 0, 

5(2n + 5- 1) 



(4) 



[z + n)s 



1 + 



2z 



+ 0(0 



|arg z| < 7T — e, e > 0, 



we deduce that 1 a ^{q; z) is an entire function in 2 and is analytic on the domain {q £ C; > 
— 1}. For b < 0, we recall that there exists 6 > such that V'(^) = and thus ipe(u) = ip(9 + u). 
In this case, by setting 6 a = — , it is shown in 20] that there exists a positive constant Ce a such 
that 

2a,v fa") ~ Ce a x 9 l a ^ g (x a ) as x -> 00. 



We also introduce the function M a ^fi(q\ x a ) denned by 

(5) KM* x °) = ^(gi * a ) ~ C eci x e T{q T ^ a) l a ^ (q + 6 a ; x a ) , K(x)>0. 

Moreover, if we assume that there exists [5 G [0, 1] and a constant > such that lim^oo ip(u) / u 1+/3 
ag, then Cg Q is characterized by 



Co 



v{i-e a ) (e a -i)\ 



if 9 a is a positive integer, 



a a 



.Shi. {k+e a )^{ak) 



1 lfc=i e k^{ak+9 a ) < otnerwise, 



where -B 7 stands for the Euler-Mascheroni constant. We recall, also from [20(], that, for r, x > 0, 



Ex 



e rJ ° 



J Q ^ (r/) - Ce^xfl^ (rx a ) . 



We deduce from Theorems 11.21 and 11.41 the following. 
Corollary 3.2. Xei q > and < x < a. Then, 



E, 



and 



1 + X?(a) 



where T/*\ = inf{u > 0; X u = a(l + x u ) a }- We also deduce that 



E, 



6> 



Moreover, assume b > and set /? = aAx and 7 = <fi(q) . Then 

1 \ 7 ^a,V>7 (a'X, 



-uT z 
3 y- 1 a 



fee/ UftxW, 



, < - < a, 



E„ 



7. (a 1 ! ■ 
a' I 1 / 



5-^-, < x < a. 



Finally, if b < and < 9 < a, we have 

E r \e-i T o 



a,if>,6 \x' ,XX 



a, 



6;x^ a 



Remark 3.3. From the strong Markov property and the absence of positive jumps, we easily get 
that first passage times above for the processes U and Z are infinitely divisible random variables. 
Hence, we obtain from Corollarv l3.2l that the functions ([3|) and ([5]) are Laplace transforms, with 
respect to the parameter q, of infinitely divisible distributions concentrated on the positive real 
line. 

We end up by investigating some special cases which allow to make some connections between 
the power series introduced and some well-known or new special functions. 

3.2.1. The confluent hypergeometric functions. We first consider a Brownian motion with drift 
—v, i.e. ip(u) = \u 2 — uu. Setting a = 2, we have 9 = 2v and therefore we assume v < 1. 
Its associated semi-stable process is well known to be a Bessel process of index v and thus 
the associated Ornstein-Uhlenbeck process is, in the case n = 2u + 1 6 N, the radial norm of 
re-dimensional Ornstein-Uhlenbeck process. We get 

T 2 ,i,{x) = (x/2y/ 2 T(-u + 1)I_„ (yfa 

where I u (x) = X^^o w ir(^+ w +i) stands for the modified Bessel function of index v, see e.g. 
5.], and 



17, 



/ x 2 

^2,ip(q;x 2 ) = $lq,i-v,— 

( x 2 

where Q(q,v, x) = Yl^Lo iy) m! ^" stands for the confluent hypergeometric function of the first 
kind, see e.g. [13, 9.9]. Using the asymptotic behavior of the Bessel function 

e x 

I u (x) ~ —== as x — > oo, 
\2ttx 

we deduce that Civ = — V y(v) • Hence, 

Kf t 2^ (*( 1 X<2 \ 2^ r (-^) r (^ + ^)^ ( i X 2 



r(i/) V 2 

2 

where A.(q,u + 1, ^-) is the confluent hypergeometric of the second kind. We mention that, in 
this case, the results of Corollary 13.21 are well-known and can be found in Matsumoto and Yor 
[l8| and in Borodin and Salminen 0, II. 8. 2]. 



3.2.2. Some generalization of the Mittag-Leffler function. Patie 2l( introduced a new parametric 
family of one-sided Levy processes which are characterized by the following Laplace exponents, 
for any 1 < a < 2, and 7 > 1 — a, 

(6) ^ 1 {u) = -{{u + 1 -l) a -{ 1 -l) a ). 

a 
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Its characteristic triplet are a = 0, v{dy) = p(2-a) y < 0, and & 7 = (7)0 (T (7 — 1 + 

a) — T(7 — 1)) where T(A) = Y(Xy ^ s ^ e digamma function. In particular, if 70 denotes the zero 
of the function 7 — ► 6 7 , then for 7 > 70 € (1 — a, 0), b > 0. 



TTie case 7 = 0.© reduces to ip(u) = ^(u — l) a . Observe that 6 = 1, if^'(l) = and 

! r(q(ra + l)) 
T{a) 

The series ([3]) can be written as follows 

T\4,{q\x) = T{a-l)M q aa _ l {ax) 



where 



stands for the Mittag-Leffler function of parameter a, (3,q > 0, which was introduced by Prab- 
hakar [23|]. Moreover, we have, see e.g. [20 ]. 

^ —, ~ —e x x l ~ f3 l{x a ) ass- > 00, 

^ o r(an + /3) a 1 ; 

with Z a slowly varying function at infinity. Thus, Ci = -§r and 

a 

As concluding remarks, we first mention that in the diffusion case, i.e. when (U,Q X ) is the 
Ornstein-Uhlenbeck process associated to a Bessel process, see 13.2. H the law of the first pas- 
sage time above can be expressed as an infinite convolution of exponential distributions with 
parametersgiven by the sequence of positive zeros of the confluent hyper geometric function, 
see Kent [131 ] for more details. Beside this case, we do not know whether such a representation 
is available. For instance, the location of the zeros of the generalized Mittag-Leffler functions, 
considered in the second example treated above, is still an open problem, see e.g. Craven and 
Csordas[8(. 
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